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Abstract: Proton and carbon ion therapy is an emerging technique used for the treatment of solid 
cancers. The monitoring of the dose delivered during such treatments is still a matter of research. 
A possible technique exploits the information provided by single photon emission from nuclear 
decays induced by the irradiation. This paper reports the measurements of the spectrum and rate 
of such photons produced from the interaction of a 80 MeV / u fully stripped carbon ion beam at 
the Laboratori Nazionali del Sud of INFN, Catania, with a Poly-methyl methacrylate target. The 
differential production rate for photons with energy E > 2 MeV and emitted at 90" is found to be 
dNy/{dNcdQ.) = (2.92±0.19) x lO^^sr-i. 
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Introduction 

In the last decade, the use of proton and carbon beams has become more and more widespread as 
an effective therapy for the treatment of solid cancer (hadrontherapy). Due to their very favorable 
profile of the release dose in tissue, the hadron beams can be very effective in destroying the tumor 
and sparing the adjacent healthy tissue in comparison to the standard X-ray based treatment OU. On 
the other hand, the space selectivity of the hadrontherapy asks for a new approach to the delivered 
dose monitoring. 

The uncertainty on the position of the dose released in hadrontherapy treatment can be due to vari- 
ous factors: calibration of the Computed Tomography (CT) images, possible morphologic changes 
occurring between CT and treatment, patient mis-positioning and organ motion during treatment 
itself. All these effects give an overall uncertainty of the order of few milUmeters that can be larger 
than the dimension of the peak of the dose release (Bragg Peak). A precise monitoring of the dose is 
then essential for a good quality control of the treatment. Furthermore, the dose monitoring would 
be particularly useful if provided during the treatment (in-beam monitoring) in order to provide a 
fast feedback to the beam. 

Several methods have been developed to determine the Bragg Peak position online by exploiting 
the secondary particle production induced by the hadron beam. Since the in^adiation of tissues 
with hadron beams produces nuclear excitations followed by photon emissions from de-excitation 
within few ns (prompt photons), both rate and production region could be used to monitor the 
dose 

In this paper we present the measurement of secondary particles produced by the Poly-methyl 
methacrylate (PMMA) during an irradiation with carbon ions. In section [I], we describe the setup 
of the experiment performed at the Laboratori Nazionali del Sud (LNS) Istituto Nazionale di Fisica 
Nucleare (INFN) di Catania with cerium-doped lutetium-yttirum ortho-siUcate (LYSO) crystals 
for the detection of prompt photons emitted in the interaction of a 80 MeV/u fully stripped carbon 
beam with a PMMA target. The obtained results on the measured spectra and rates are then reported 
in sections and Q|. 
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1. Experimental setup 



The experimental setup is sliown in Fig. A 4 x 4 x 4 cm^ PMMA target is placed on 960 MeV 
(80 MeV/u) fully stripped ion beam. The beam rate, ranging from hundred kHz to ~ 2MHz, is 
monitored with a 1 . 1 mm thin scintillator on the beam line read with two photomultiplier tubes 
(PMT) (Hamamatsu 10583) put in coincidence (Start Counter) and placed at 17 cm from the 
PMMA. 



Beam 




Figure 1. Schematic view of the experimental setup; the acquisition is triggered by the Start Counter in 
coincidence with the LYSO. In the picture of the detectors is also shown. 




Prompt photons are discriminated by measuring a time difference between the beam impact 
and the photon detection consistent with a particle emitted istantaneously and traveling at the speed 
of light. The main background to prompt photons is due to neutrons, that are also emitted by 
nuclear interactions of the carbon beam in the PMMA. Since such neutrons are non relativistic, 
their arrival time is a good discriminator between photons and neutrons, hence the time resolution 
becomes critical in the signal selection. Therefore, we have chosen the LYSO crystals and tested 
their performance in discriminating between prompt photons and background neutrons. LYSO 
crystals have high light output, high density, and fast decay time (as an indication see Tab. HJ) and 
are commonly exploited in PET imaging, although in smaller sizes. 

An array of 4 LYSO crystals, each measuring 1.5 x 1.5 x 12 cm^, is placed at 90° with respect 
to the beam line, at 74 cm from the PMMA center. The scintillation light of the crystals is detected 
with a PMT (EMI 9814B) triggered in coincidence with the Start Counter. A 12-bit QDC (Caen 
792N) and a 19-bit TDC (Caen 1 190B) provide the measurements of both the particles' energy and 
arrival time. 

1.1 Calibration and Simulation 

In order to calibrate the LYSO detector we used ^^Na and '^"Co sources: fitting gaussian peaks 
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LYSO Characteristics 




Eiiective Atomic Number 


66 


Density {g/crsv') 


1 A 


Radiation Length (cm) 


LIO 


Decay Constant (ns) 


40-44 


Peak Emission (nm) 


428 


Light Yield % Nal (Tl) 


75 


Index of Refraction 


L82 



Table 1. Properties of LYSO crystals. 



for the observed energy spectra, linearity is verified up to the energy of the second Co peak 
(1.33 MeV). At higher energies, when the pair-production mechanism becomes relevant in the 
shower development inside the crystals, there is a non-negligible probability to escape the detector 
for one or both 511 keV photons produced by the positron annihilation, mechanism known as 
single- and double- escape respectively. When such phenomena occur, the output signal appears as 
a superposition of three gaussian structures and, therefore, the measured energy is biased and its 
resolution degrades. 

This effect was simulated using the GATE framework [0], a tool dedicated to medical imaging, 
radiotherapy and hadrontherapy, based on the GEANT4 Monte Carlo (MC) code Hl.EIl- 




true E (MaV) true E IMeV] 

Figure 2. Results obtained using GATE simulation. Left (right): mean values (resolution/energy) of the 
energy spectra recorded in LYSO detector, as a function of the photon energy. The red line in the left picture 
represents the bisector line. 

We simulated the response of the experimental setup to monochromatic photons of several 
energies. The simulated distribution of the collected photoelectrons have been fitted with gaussian 
functions with exponential tails, and the mean values for energies below 1 MeV were used for 
energy calibration, thus following the same procedure used for data analysis. Using simulated pho- 
tons of energies above 1 MeV we estimate the bias on the energy measurement and the degradation 
of the resolution due to single- and double- escapes. The resolution obtained for the calibration 
peaks determines the number of photons to be used in the MC simulation. The results are shown 



-3- 



in Fig. 0. The left plot shows the mean value of the energy reconstructed simulating the response 
of the experimental setup as a function of the true energy; the comparison with the bisector (red 
line) shows the ~ 10% bias on the energy measurement due to single- and double- escapes for 
E > I MeV. On the right plot, the same results are shown for the a/E parameter as a function of 
the true energy. A deterioration of the resolution for £ > 1 MeV is observed. 

1.2 Event timing 

In order to select signal events, the time difference between the energy deposition in the LYSO 
detector (Tlyso) and arrival time of the ^^C ions on the Start Counter (Tsc), referred to as At, is 
evaluated. When both photo-multiphers of the Start Counter yield a signal, the smallest absolute 
value of At is chosen. 

The number of carbon ions reaching the target and the electronics inefficiency due to the 
acquisition dead-time are needed to estimate the prompt photon rate. The number of carbon ions 
reaching the PMMA in a given time interval is computed by counting the number of signals given 
by the Start Counter (Nsc) within randomly-triggered time-windows of = 2 \xs. Measuring the 
total run time (Ttot) and the number of time windows considered (Nw), the number of carbon ions 
is estimated as 

£sc N„ * T„ 

The Start Counter efficiency esc = (96 ±1)% has been estimated by exploiting the two-sided 
PMT readout with negligible dark counts: esc = £sc\ • £sc2 = ^1^2/ (^1 ■^2)> with N\^2 are the single 
PMT counts and their coincidences. 

To con^ect Nc for the dead-time inefficiency, eoT has been estimated from the total acquisition 
time {Tdead) as: 

eDT = l-^. (1.2) 

The measured values of e^T at an average carbon ion rate Rc = Nc/Ttot ranges from 70% at 
0.6 MHz up to 47% at 2 MHz. This efficiency correction was then applied to data. 

2. Prompt photon energy spectrum measurements 

The con^elation between the reconstructed energy (E) and the measured At is shown in Fig. ^ (Left). 
Three populations are present: a horizontal band from LYSO intrinsic noise; an almost vertical 
band due to the signal from the prompt photons; a diffused cloud mainly due to neutrons at At 
values larger than those of the prompt photons population. Events are shown in a ±23 ns window 
around the average At values for prompt photons to avoid events associated with the previous and 
subsequent cyclotron cycles. The measured horizontal offset in the At distribution for prompt 
photons is due to the electronic setup. A secondary component, few ns before the most populated 
band, is visible: this distribution is due to the prompt photons produced inside the Start Counter. 
The earlier arrival time of these events is due to their shorter travel. The shape of the prompt 
photon band is instead due to the time slewing effect induced by the front-end electronics fixed 
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Figure 3. Left: Calibrated energy released in LYSO crystal as a function of the arrival time, Af. Three 
component are present; an horizontal band due to the LYSO intrinsic noise; an almost vertical band due to 
the signal from the prompt photons; a diffused cloud mainly due to neutrons at values larger than those 
of the prompt photons population. Right: estimated time slewing correction. 



voltage threshold. This effect can be corrected for by fitting the A/ distribution in bins of E, as 
shown in Fig. ^ (right), thus obtaining the correction C{E). 




AT Corrected [ns] I ^"""^'Z I""^! 



Figure 4. Left: calibrated energy released in LYSO crystal as a function of the arrival time corrected taking 
into account the slewing effect. Right: measured time resolution as a function of the calibrated energy. 

Fig. (Left) shows the energy spectrum as a function of the corrected time (A/^orr = A/ — 
C{E)). Fig. fright shows the time resolution {Gai) as a function of the energy. A resolution better 
than 300 ps is achieved for energies greater than 3 MeV. 

In order to measure the prompt photon energy spectrum originating from the PMMA, the num- 
ber of prompt photons in each energy bin is extracted from the fit to the distribution of Atcorr/ f^Ar^rr 
(pull) in that bin: the fit uses a polinomial function and a gaussian distribution, fixing the width at 
one and the mean at zero. Fig. shows an example of the pull distribution obtained for one energy 
bin. The area under the gaussian provides the number of prompt photons in a given energy bin. 
The resulting measured spectrum, normalized to the number of incident carbon ions and corrected 
for Edt is shown in Fig. ^ (left). 
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Figure 5. Example of a pull distribution and the corresponding fit to extract the number of prompt photons. 



The experiment setup has been simulated with GATE to obtain the expected energy spectrum : 
we used G4QMDReaction to model ions' inelastic processes and the interactions between '^C and 
the PMMA. The simulated spectrum of the prompt photon emitted in the geometrical acceptance 
of the LYSO crystal is shown in Fig. § (Right). Without the detector effects, two components are 
visible: a continuum component and few excitation peaks. The two most prominent structures at 
E = 4.44 MeV and 6.13 MeV come from ^^C* and ^^O* excited states decays respectively. 




Energy [MeV] Energy [MeV] 



Figure 6. Left: Measured prompt photons spectra obtained with the LYSO detector The fit of the peak gives 
a fraction of fuf^ = (13.9±0.6)% for prompt photons over carbon ions at 4.44 MeV. Right: Monte Carlo 
energy spectrum of the prompt photons in the detector acceptance without detector response simulation. 



Fitting the data using exponential function for the continuum contribution and a gaussian dis- 
tribution for the main ^^C* peak we obtain for this peak a mean value and a resolution compatible 
within errors with the expectations from the simulation: Ge ~ 0.5 MeV, for E = 4.44 MeV (see 
Fig. From this fit we measure a contribution from the ^^C* decays Une to the overall rate of 
prompt photons with energies E > 2MeV equal to /n^ = (13.9 ±0.6)%. 

The MC true energy distribution has been folded with the detector response (detector and Start 
Counter efficiency, resolution, and bias), using the simulations described in Sec. [I]. The folded 



-6- 



MC spectrum is compared with the measured one in Fig. [7|, with the spectra normalized to the 
number of incident carbons. The comparison shows the level of agreement between data and 
MonteCarlo simulation: both the normalization and the fraction of E = 4.44 MeV photons require 
further investigation. 
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Figure 7. Data-MC comparison of the observed energy spectrum of the prompt photons. All spectra are 
normalized to the number of incident carbon ions. 




3. Rate measurements 

The ultimate goal of the studies on prompt photons is the measurement of the number of carbon ions 
delivered to the patient. To understand the potentiality of LYSO detectors as dosimeters for prompt 
photons measurements, we select only events with E > 2 MeV rejecting the background from 
natural radioactivity of the LYSo crystal. For such events we proceed with a side-band subtraction 
on the prompt photon At^orr distribution and extract the number of prompt photons, Nmeas- For runs 
with different carbon ion rates we then measure the ratio between the rates of observed prompt 
photons (Rprompt) and carbon ions (Re): 

p Rprompt Nmeas ^ 

prompt SdT^ScNc 

with Edt, £sc^ and Nc the dead-time efficiency, the Start Counter efficiency, and the number of 
carbon ions respectively, measured on data (Sec. [I]). Fig. ^ shows that Fprompt does not dependent 
on the carbon rate and averages to 

Fprompt = (3.04 ±0.01,,,,, ±0.20,,,) X 10-^ (3.2) 

The systematic error was introduced to account for the dispersion of the values which is well 
above the statistical fluctuations. Under the hypothesis of a flat distribution, it was estimated as the 
semi-dipersion divided by \/3. 

We then measured the differential prodution rate: 
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Figure 8. Rate of prompt photons with E > 2 MeV normahzed to the rate of carbon ions as a function of 
the latter. The red fit line is shown with the error band (both statistical and systematic errors are taken). 



{E > 2MeV, 6=90")= ^""""p' (3.3) 

dNcdO. Smeas^LYSO 

with D.LYSO = 1 -3 X 10~^ the solid angle covered by the detector and Emeas = (8 1 .3 ± 2.5) % the 
detection efficiency estimated with the Monte Carlo simulation (Sec. [I]). To evaluate the uncertainty 
due to the assumed energy distribution of the photons, this efficiency has been computed both using 
a flat energy distribution and using the simulated MC spectrum. The semi-difference between the 
two results is then assumed as the uncertainty. 

The final result for the differential rate of prompt photons is: 

dNcdO. 
Conclusions 



[E > 2MeV,0 =90'0 = (2.92±0.01,,«, ±0.1%) x 10"V"\ (3.4) 



The measurement of the energy spectrum and rate of prompt photons produced in the interaction of 
80 MeV/u carbon ions with a PMMA target has been reported. These results were obtained with 
a LYSO detector with a time resolution of 300 ps for photons with energies above 3 MeV. The 
excellent time resolution was used to discriminate the neutrons and to determine a precise energy 
spectrum of the prompt photons. In particular, we measure the fraction of prompt photons with 
energies E>1 MeV coming from the E = AAA MeV line to be fn^ = (13.9 ±0.6)%. 

We also measure the rate of photons per carbon ions triggered to be Fprompt = (3.04±0.0l5(flf ± 
O.lOsys) X 10^^. We also proved Fprompt is largely independent, after dead-time corrections, from 
the carbon rate in the range between 0.6 and 2 MHz explored in this experiment. 
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Finally, we measured the differential production rate to be dNy/dNcdQ. {E > 2 MeV, 6 = 
90 °) = (2.92±0.19) X 10-2 sr'^ 
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